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This positron emission tomography (PET) study aimed to further define selectivity of [ 11 C]Ro15-4513 
binding to the GABARa5 relative to the GABARal benzodiazepine receptor subtype. The impact of 
Zolpidem, a GABARal -selective agonist, on [ 11 C]Ro15-4513, which shows selectivity for GABARa5, 
and the nonselective benzodiazepine ligand [ 11 C]flumazenil binding was assessed in humans. 
Compartmental modelling of the kinetics of [ 11 C]Ro15-4513 time-activity curves was used to 
describe distribution volume (V T ) differences in regions populated by different GABA receptor 
subtypes. Those with low a5 were best fitted by one-tissue compartment models; and those with 
high a5 required a more complex model. The heterogeneity between brain regions suggested 
spectral analysis as a more appropriate method to quantify binding as it does not a priori specify 
compartments. Spectral analysis revealed that Zolpidem caused a significant V T decrease (~ 10%) in 
[ 11 C]flumazenil, but no decrease in [ 11 C]Ro15-4513 binding. Further analysis of [ 11 C]Ro15-4513 
kinetics revealed additional frequency components present in regions containing both art and «5 
subtypes compared with those containing only ort. Zolpidem reduced one component (mean ±s.d.: 
71% ±41%), presumed to reflect a1-subtype binding, but not another (13% ±22%), presumed to 
reflect ot5. The proposed method for [ 11 C]Ro15-4513 analysis may allow more accurate selective 
binding assays and estimation of drug occupancy for other nonselective ligands. 
Journal of Cerebral Blood Flow & Metabolism (2012) 32, 731-744; doi:10.1038/jcbfm. 2011. 177; published online 
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Introduction 

Imaging studies of benzodiazepine receptors over the 
past 20 years have revealed different levels of 
binding in patients with anxiety and panic disorders, 
alcoholism, schizophrenia, epilepsy, and other dis- 
orders (Malizia et al, 1998; Lingford-Hughes et al, 
2010; Busatto et al, 1997; Hammers et al, 2001). 
There is increasing knowledge about the six different 
subtypes of the GABA-benzodiazepine receptor, their 
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locations in the brain, and their functional implica- 
tions. For instance, knock-in and knockout genetic 
manipulations in mice have shown that the GABA- 
benzodiazepine receptor subtype containing the al 
subtype is important in sedation and sleep, the a2/3 
subtype is involved in anxiety, and the <x5 in learning 
and memory (Mohler et al, 2002). In the light of 
these discoveries, it is important to study the role of 
these receptor subtypes in humans. Neuroimaging, 
in particular positron emission tomography (PET) 
and single-photon emission tomography (SPET) 
imaging with receptor-selective ligands, remains 
the only way of quantifying benzodiazepine recep- 
tors in the human brain in vivo. 

[ 11 C]flumazenil with PET is the most commonly 
used in vivo marker of benzodiazepine receptor 
binding and its iodinated analogue, [ 123 I]iomazenil, 
can be used with SPET. Flumazenil is an antagonist 
that binds with high affinity to benzodiazepine 
receptors of the al, a2, a3, or a5 (Ki ~lnmol/L) 
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subtypes and with lower affinity to the a4 or a6 
subtypes (Ki ~150nmol/L (Sieghart, 1995)). Until 
recently, there was no way of preferentially visualis- 
ing particular subtypes. Rol5-4513 shows 10- to 15- 
fold selectivity for the a5 subtype over others in vitro 
(Ki 0.07nmol/L versus 1 to 5nmol/L), and in vivo 
[ 11 C]Rol5-4513 PET provides a relatively selective 
marker for the GABARa5 subtype (Hadingham et al, 
1993; Lingford-Hughes et al, 2002; Maeda et al, 
2003). These studies, among others, show that in vivo 
[ 11 C]Rol5-4513 uptake targets the limbic system, 
with the anterior cingulate cortex, ventral striatum, 
and hippocampus showing particularly high levels 
(Onoe et al, 1996), whereas high [ 11 C]flumazenil 
binding is seen in all cortical areas, with highest 
levels in the occipital cortex. Therefore, the limbic 
pattern of [ 11 C]Rol5-4513 binding is likely to be 
relevant to learning, memory, reward, and addiction. 
For example, we have shown that an inverse agonist 
at the a5 subtype can reverse the memory-impairing 
effects of acute alcohol consumption (Nutt et al, 
2007), and ["C]Rol5-4513 binding is related to 
performance on a verbal memory task in alcohol 
dependence (Lingford-Hughes et al, 2010). The a5 
subtype has also been implicated in the reinforc- 
ing effects of alcohol in mice (Stephens et al, 2005), 
and we have shown a reduction in [ 1:1 C]Rol5-4513 
binding in the nucleus accumbens in alcohol depen- 
dence compared with healthy controls (Lingford- 
Hughes et al, 2010). 

Although the distribution of ["C]Rol5-4513 bind- 
ing is consistent with relatively greater labelling of 
the a5 subtype, the affinity of [ 11 C]Rol5-4513 is only 
an order of magnitude greater for the a5 over other 
subtypes. Given the abundance of the al subtype 
throughout the brain, it is important to assess the 
contribution of al binding to ["C]Rol5-4513 bind- 
ing. This study compares the binding of [ 11 C]Rol5- 
4513 with that of ["CJflumazenil in the presence and 
absence of Zolpidem, a compound with selectivity 
for the al receptor subtype that is > 5,000-fold over 
that for the a5 subtype (Hanson et al, 2008; Langer 
et al, 1992). Unlike classical benzodiazepines, the 
sedative/hypnotic effect of Zolpidem occurs at 
much lower doses than the other pharmacological 
effects attributed to benzodiazepine-site action, such 
as anticonvulsant activity and muscle relaxation 
(Sanger, 2004). 

After predosing with Zolpidem or placebo, 
one group of healthy volunteers was scanned with 
["0^015-4513 and another group with ["CJfluma- 
zenil. Our hypothesis was that, because [ 11 C]fluma- 
zenil labels the al subtype, its binding would 
be reduced by Zolpidem. In contrast, we hypo- 
thesised that ["C]Rol5-4513 binding would not be 
altered by Zolpidem in a5-rich regions such as 
the hippocampus, whereas in al-rich regions with 
limited or no a5, such as the cerebral cortex 
and thalamus, ["C1R015-4513 binding might be 
reduced. To further resolve contributions of al and 
a5 subtypes to the [ 11 C]Rol5-4513 signal, we applied 



spectral analysis (Cunningham and Jones, 1993) to 
both [ 11 C]Rol5-4513 and ["Qflumazenil data sets to 
deduce the spectral band of the GABA-benzodiaze- 
pine al subtype. The range obtained was used 
to determine the occupancy by Zolpidem of the 
GAB A al subtype from the ["C]Rol5-4513 data using 
the Lassen plot (Lassen et al, 1995; Cunningham 
et al, 2010). 

Materials and methods 

This study was approved by a NHS Research Ethics 
Committee, the Administration of Radioactive Substances 
Advisory Committee and local NHS Research and Devel- 
opment. It was conducted in accordance with Good 
Clinical Practice — International Conference on Harmonisa- 
tion guidelines. 

Design 

There were two separate randomised double-blind placebo- 
controlled studies, each with six healthy volunteer parti- 
cipants. The study procedures were identical apart from 
the PET ligand used: [ 11 C]flumazenil for one study and 
[ 11 C]Rol5-4513 for the other. Positron emission tomo- 
graphy scans were performed twice for each participant, 
once after Zolpidem and once after placebo, with at least a 
week between scans. Binding of the ligand was compared 
between the two conditions. 

Subjects 

Healthy male volunteers (age mean + s.d.: [ 11 C]Rol5-4513: 
44 ± 6 years; [ 11 C]flumazenil: 43 ± 4 years) underwent 2 
PET scans at least a week apart, after Zolpidem or placebo, 
the order of which was randomised. They were all 
screened with urinalysis for illicit drugs of abuse, had no 
current or previous significant physical or mental dis- 
orders, had taken no psychotropic drugs in the last month, 
nor were regular benzodiazepine users. 

Zolpidem and Measures 

Zolpidem (20 mg) or an identical placebo was administered 
90 minutes before radioligand injection and at the start of 
the PET scan. Plasma Zolpidem levels were measured at the 
time of PET ligand injection. 

The following tasks provided a measure of the effect 
of Zolpidem on brain function. Saccadic eye movements 
to a target were measured using electrooculography as 
described previously (Lingford-Hughes et al, 2005). The 
subject was asked to fixate on a red light as it appeared 
on the screen mounted 67 cm in front of him. Data were 
then collected for 48 eye movements of 15° to 40° and a 
main sequence curve obtained, the measurement for peak 
velocity being the interpolation at 35°. Verbal memory was 
tested with learning and delayed recall of a list of 20 words 
(Nutt et al, 2007), with the list presented at the start of the 
scan and recall 30 minutes later. 
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Positron Emission Tomography Method 

["C]Rol5-4513 was synthesiser! by N-methylation of the 
JV-desmethyl derivative with [ 11 C]iodomethane. Reverse- 
phase high-performance liquid chromatography was used 
to purify the compound (Phenomenex Ultra-carb 7 ODS, 
250 x 10mm 2 ). [ 11 C]flumazenil was synthesised by methyla- 
tion of the nor-derivative with "CHJ. 

A bolus injection of [ 11 C]flumazenil (mean + s.d.: 
476 + llMBq in 1.0±0.16mL) or ["C1R015-4513 
(meanls.d.: 495±18MBq in 3.3 + 0.49 mL) was adminis- 
tered through an intravenous cannula sited in the domi- 
nant antecubital fossa vein. No significant differences in 
cold mass injected, with paired f-tests, between placebo 
and Zolpidem scans were found using either flumazenil 
(meanls.d.: 1.37 ± 0.52 /tg to 1.09 ± 0.28 /jg, f=1.28, 
P>0.25) or ["C1R015-4513 (meanls.d.: 3.22 ± 0.86 ng to 
4.30 + 3.71 ng, f=0.579, P>0.5). 

We used a Siemens ECAT EXACT HR+ (CTI/Siemens, 
model 962; Knoxville, TN, USA) scanner with an axial 
field of view of 15.5 cm. Subjects were positioned such that 
the transaxial images were parallel to the intercommissural 
line. A total of 63 transaxial image planes were acquired 
as 2.42-mm slices with a reconstructed axial resolution 
of 5.4mm and a transaxial resolution of 5.6 mm (Spinks 
et al, 2000). A 10-minute transmission scan was acquired 
before every emission using a single rotating photon 
point source of 137 Cs for subsequent attenuation and scatter 
correction. 

[ 11 C]Rol5-4513 scans comprised 24 dynamic time frames 
(1 x 30, 4 x 15, 4 x 60, 2 x 150, 10 x 300, 3 x 600 seconds) 
and [ 11 C]flumazenil scans comprised 20 time frames (1 x 30, 
4 x 15, 4 x 60, 2 x 150, 2 x 300, 7 x 600 seconds) of data 
which, after corrections for attenuation, random coinci- 
dences, and scatter, were all reconstructed using filtered 
backprojection. The three-dimensional emission data were 
Fourier-rebinned (Defrise et al, 1997) into two-dimensional 
data sets to increase the speed of reconstruction. 

Each subject had a radial arterial cannula inserted in the 
nondominant wrist to allow continuous counting of blood 
radioactivity concentration for the first 15 minutes of the 
experiment using a bismuth germanate counter. Discrete 
samples were also taken 4, 6, 8, 10, 20, 35, 50, 65, 80, and 
90 minutes after injection. An aliquot of each discrete 
sample was rapidly centrifuged to obtain corresponding 
plasma and radioactivity concentrations. The continuous 
blood counts were corrected using the estimated plasma/ 
blood ratio to derive the plasma radioactivity profile. The 
plasma was further analysed for radiolabeled metabolites 
using a semiautomated analysis system, consisting of an 
on-line solid-phase sample enrichment column and high- 
performance liquid chromatography. The amount of 
[ 11 C]Rol5-4513 or [ 11 C]flumazenil in the plasma at a given 
time point was calculated as a percentage of all radioactive 
components separated by high-performance liquid chro- 
matography added to the radioactivity present in the solid- 
phase extraction fraction, to produce a final plasma input 
function derived for radioligand-binding quantification. 
These calculations were carried out using Clickfit, in- 
house software running in Matlab (The Mathworks, Natick, 
MA, USA). 



All subjects had Tl-weighted magnetic resonance imag- 
ing with a Philips 1.5-T Gyroscan Intera scanner (Philips, 
Best, The Netherlands) to produce structural images for 
reference. 

Image and Blood Analysis 

Reconstructed PET images were analysed on a Sun 
SPARC workstation (Sun Microsystems, Mountain View, 
CA, USA) using Analyze AVW version 8.1 (Biomedical 
Imaging Resource, Mayo Clinic, Rochester, MN, USA), 
Matlab 6 and 6.5 (The Mathworks), and SPM5 (available 
through http://www.fil.ion.ucl.ac.uk/spm/; FIL, 12 Queen 
Square, London, UK). 

Regions of interest were placed on the images using a 
probabilistic map which defined 83 regions (Hammers 
et al, 2007). A list of these regions is given in Appendix 1. 
Each structural magnetic resonance image was coregistered 
to the PET add image, a sum of all counts acquired during 
the scan, using a rigid-body technique. This coregistered 
magnetic resonance imaging was then normalised to the 
'Montreal Neurological Institute' space (International 
Consortium for Brain Mapping (ICBM)/MNI) using a 
bias-corrected segmentation in SPM5. Using the reverse 
parameters of this normalisation, the probabilistic map 
could be warped to fit the coregistered structural magnetic 
resonance imaging and hence the add image and the original 
PET dynamic. The goodness of fit of each object map to the 
individual brain was checked visually before regional 
sampling using the region of interest tool in Analyze 8.1. 
From the probabilistic map, we identified regions with 
wide differences in al and a5 distribution: the cerebellum, 
occipital cortex, and thalamus, which contain high al levels 
and very low a5 levels, and limbic regions such as the 
hippocampus, nucleus accumbens, amygdala, and anterior 
cingulate cortex with high a5 levels and lower al levels 
(Fritschy and Mohler, 1995; Lingford-Hughes et al, 2002). 

Modelling 

Parametric maps of total volume of distribution [V T ) were 
created for both [ 11 C]Rol5-4513 and ["C] flumazenil images 
using receptor parametric mapping software written in 
Matlab (RPM, MRCCU (Cunningham and Jones, 1993)). 
Values for mean V T per voxel in each region of interest and 
in combined regions were then sampled in Analyze 8.1, to 
compare in vivo binding with and without Zolpidem, for 
the two different PET ligands. 

The regional kinetics of uptake of [ 11 C]Rol5-4513 only 
were computed from brain and plasma time-activity curves 
(TACs) using one-tissue and two-tissue compartment 
models, using Clickfit. Having estimated the delay, the 
models used fixed blood volume (V b ) of 0.06, calculated 
after a preliminary analysis of variable V b and averaged 
over 83 regions of interest. Goodness of fit for each model 
in every region was represented using Akaike information 
criteria. V T of [ 11 C]Rol5-4513 was then calculated from 
the estimation of kinetic parameters in the two models 
using the following equations: one-tissue compartment 
V T = KJk 2 ; two-tissue compartment V T = || f 1 + |M . 
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Spectral Analysis 

To further elucidate the pharmacokinetics of [ 11 C]Rol5- 
4513, Clickfit was again used to carry out spectral analysis 
of each regional TAC collected from the original PET 
dynamics. Spectral analysis (Cunningham and Jones, 1993; 
Turkheimer et al, 1994) convolves the arterially derived 
metabolite-corrected plasma input function with a first- 
order poly-exponential, without specifying the number of 
tissue compartments, to fit the TAC for any region, under 
the assumptions of system linearity and time invariance. 
The impulse response function (IRF) is described as 
IRF(t) = Yli 8/6^*1 where a and [J are real, nonnegative 
numbers. Overall, 100 discrete values of /J on a logarithmic 
scale (Turkheimer et al, 1994), corresponding to 100 
different basis functions, were chosen, between 0.000631 
to 1.0 seconds -1 , a lower limit shown to show the greatest 
interclass correlation coefficient (Barros et al, 2010). The 
'spectrum' of peaks with height i, can be used to infer 
pharmacological characteristics of the ligand and its 
binding properties. Given that the area under the IRF 
function is the total volume of distribution, the position 
and height of spectral peaks can be used to estimate V T 
using V T = Ylk=i i a k/Pk ~ fy, where n is the total number of 
peaks, a k the Ath peak height, p k is kth peak position, and X 
the decay constant. 

In the case of two specific binding sites, the total volume 
of distribution V T can be written as V T = V S1 + V S2 + Vnd, the 
sum of the specific distribution volumes, V S1 and V S2 , and 
the nondisplaceable distribution volume V ND . Particularly, 
for [ 11 C]Rol5-4513, the V S1 relative to the specific binding 
to GABARot5 and V S2 relative to the specific binding to 
GABARal may be determined comparing the contribution 
of different kinetics to the IRF in a5-rich regions with those 
in region with little or no or.5 density, before and after 
blocking of the al-binding signal by Zolpidem. As the IRF 
is linear, V s may be estimated for particular regions of the 
spectrum, Vs = Yl'k=i ( a k/Pk ~ where 1 and /are integers, 
j>i, representing /^-values defined by pharmacological and 
regional inference. The ranges of j8 chosen to estimate Vs 
for different subtypes are described in detail in the 
'Results' section. All peaks with [i> 0.001 seconds -1 were, 
attributed to their fast curve shape, assumed to represent 
blood volume and nonspecific binding, and were therefore 
excluded with a low-pass filter. A similar band-pass 
approach to spectral analysis, based on regionally in- 
formed frequency allocation, has been used to investigate 
cerebral protein synthesis using L-[l- 11 C]leucine PET by 
Veronese et al, (2010). 

Occupancy Plots 

Occupancy data were derived graphically, as described by 
Lassen et al (1995) (Cunningham et al, 2010), using the 
vectors of the regional volumes of distribution in baseline, 
V$, and postdrug, Vg, as follows: V b T -V$=O d (WVmJ, 
assuming that the V ND and O d , the receptor occupancy of 
the blocking drug, were constant across brain regions. 

If the imaging ligand has two binding sites and the 
blocking drug is selective for only one of these binding 
sites, V* and V d may be written as follows: 



v b T =v b sl +v b S2 + v m 

Vt =K + V S2 + Vnd = + V| 2 (l - O d ) + V ND 

Therefore, we can draw two plots for the fractional 
volumes of distribution defined above: 

v| a -v| 1 = o 
v| 2 -v| 2 = o d v| 2 

The specific volume of distribution of [ 11 C]Rol5-4513 
relative to each subtype was derived using spectral analysis 
with the assumption that a given spectral range corre- 
sponded to binding to a single site. In addition, we made the 
further assumption that, although the blocking drug will 
change the positions of the spectral peaks, the volume of 
distribution of the imaging ligand relative to the binding to 
the site not targeted by the cold drug will be the same in the 
baseline and postdrug PET studies. Both assumptions were 
tested using Zolpidem blocking data, as reported in the 
'Results' section. The 83 region atlas (Appendix 1) provided 
a large number of regional volumes of distributions for use 
in the plots. The gradient of a linear regression through the 
points on the axes, x=(V|), y=(V|— Vs), estimated the 
Zolpidem occupancy O d at the given receptor subtype. 

Statistical Analysis 

Two-tailed Student's f-statistic was used to interrogate 
demographic variables, Zolpidem levels, as well as objec- 
tive and subjective measures between the two scanning 
groups and within groups in the presence and absence of 
Zolpidem. Mean V T estimates for the whole brain with 
and without Zolpidem were also compared with Student's 
t-test, whereas the effect of Zolpidem on binding in 
selected regions (such as the whole brain, occipital cortex, 
cerebellum, hippocampus, nucleus accumbens, frontal 
lobe, anterior cingulate, amygdala), was interrogated with 
two-way analysis of variance (ANOVA). Two-way repeated- 
measures ANOVA with P-plot graphical test and the 
Hochberg multiple comparison correction (Turkheimer 
et al, 2001) across all 83 regions was used to compare 
[ 11 C]Rol5-4513 and ["Cjflumazenil V T , with and without 
Zolpidem. One-way ANOVA was used to compare the 
estimated parameters of the different compartmental 
models, with pairwise comparisons conducted using 
Tukey's post hoc tests. To test the significance of the line 
through the regional V T estimates in occupancy plots, 
linear regressions were used. All statistical analyses were 
conducted using GraphPad Prism (version 4.02 for Win- 
dows, GraphPad Software, San Diego, CA, USA, http:// 
www.graphpad.com). 

Results 

One [ 11 C]Rol5-4513 pair of scans was not analysable 
after a subject dropped out of the study, so we report 
here the remaining five. There was no significant 
difference in Zolpidem levels in the two groups, 
although only three blood samples were successfully 
acquired in the [ 11 C]Rol5-4513 cohort (mean±s.d.: 
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["C1R015-4513 u = 3: 107.9 ± 33.1 ng/mL; ["C]fluma- 
zenil u = 6: 84.9 ±41.1 ng/mL; t=0.84; P>0.40). 
These levels are consistent with those reported in 
the literature (Salva and Costa, 1995). Zolpidem 
resulted in impaired saccadic eye movements: peak 
velocity at 90 minutes postingestion was 11.9% after 
placebo and 38.1% after Zolpidem (f=3.75; 
P< 0.005); the percentage of failed saccades was 
20.3% after placebo and 59.4% after Zolpidem 
U = 4.60; P< 0.0005). Zolpidem also impaired learn- 
ing and delayed recall of the 20-word list similarly in 
both groups. In the ["0^015-4513 group, the 
number of words recalled reduced from 6.8 ±4.2 to 
3.4 ±4.0 and in the [ 11 C]flumazenil group from 
9.2 ±5.1 to 2.2 ±1.7. 



Compartmental Modelling of ["C]Rol5-4513 

Compartmental modelling was applied to all 
[ 11 C]Rol5-4513 regional data, estimating kinetic 
parameters and goodness of fit, represented by the 
Akaike information criteria. All data for the placebo 
group are shown in Table 1, presented as mean ± s.d. 
for both models. All rate constants are expressed as 
seconds^ 1 . 

Using a one-tissue compartment model with fixed V b 
(2kbvf), estimates of Kl and k2 were consistently small 
and positive, through all regions and all individuals 
[Kl = 0.0050 ±0.0011; ic2 = 0.0015 ± 0.0008). V b was 
fixed at 0.06, calculated as the average of estimated 
V b in all 83 regions. 

The two-tissue compartment model (4kbvf), yielded, 
on average, higher estimates for Kl and k2 than did the 
one-tissue compartment model, but very low values 
for both k3 and k4, expressed as seconds 1 (4kbvf 
results: Kl = 0.0053 ± 0.0001; ic2 = 0.0019 ± 0.0001; 
k3 = 0.00042 ±0.00443; k4 = 0.00050 ± 0.00731). These 
differences were statistically significant applying a 
one-way ANOVA for each model parameter. Mean 
estimations of Kl and k2 from the four different models 



Table 1 Comparison of the kinetic rate constants obtained 
using different compartmental models to determine Rol5-4513 
binding in the brain 



Kl 



k2 



k3 



k4 



AIC 



All regions 

2kbvf 5.0 + 1.1 1.5 + 0.8 — — -89.09 + 27.8 

4kbvf 5.4 + 1.1 1.9 + 1.2 0.42 + 4.43 0.50 + 7.31 -100.76 + 28.8 



Cerebellum 

2kbvf 6.0 + 0.7 3.1 + 0.4 
4kbvf NA NA 



NA 



NA 



Hippocampus 

2kbvf 4.3 + 0.4 0.77 + 0.13 — — 
4kbvf 5.2 + 0.5 1.6 + 0.5 0.75 + 0.56 0.48 + 0.81 



-126.43 + 18.1 
NA 



-63.46 + 13.6 
-93.79 + 15.6 



AIC, Akaike information criteria; NA, not available. 

All constants are expressed as mean + s.d. seccondsx 10~ 3 . Binding in the 

cerebellum was not successfully evaluated by a two-tissue compartment model. 



were significantly different [Kl: F = 13.71; P< 0.0005, 
k2: F = 12.49; P< 0.0005), and Tukey's post hoc tests 
showed significance for both Kl and k2 [P< 0.01) in all 
comparisons of one- to two-tissue models. 

Kinetic modelling of TACs from regions rich in 
GABARal receptors gave differing results compared 
with regions rich in GABARa5 receptors. In the 
GABARal-rich cerebellum, the one-tissue compart- 
ment model fitted TACs almost as well as did the 
two-tissue compartment model (Akaike information 
criteria better in 3 of 5 individuals, average 6.8% 
more negative), but the two-tissue model produced 
consistently very small negative estimates of either 
or both k3 and jc4. However, in the GABARa5-rich 
hippocampus, there was a significantly better fit 
using the two-tissue model compared with the one- 
tissue model (Akaike information criteria better in all 
individuals, average 46.9% more negative) and 
estimates of BP (calculated using kjkj that were 
positive and small. 

Although the estimates of rate constants differed 
significantly, as did the fits of the curves, similar 
overall VrS were obtained with each model, and both 
models showed the same pattern of variability 
between individuals. 



Spectral Analysis 

[ 11 C]Rol5-4513: As previously reported, the highest 
[ 11 C]Rol5-4513 VtS were seen in limbic regions such 
as the nucleus accumbens and hippocampus (see 
Figure 1A; Lingford-Hughes et al (2002, 2010)). All 
results are presented as mean ± s.d. There was no 
reduction in ["C]Rol5-4513 V T in the presence of 
Zolpidem in the whole brain (postplacebo: 
3.69 ±0.22; postzolpidem: 3.81 ±0.47) nor among 
the combined regions (F = 0.0107; P>0.90; see 
Figure IB). Sampling all of the 83 regions from the 
probabilistic map, a repeated-measures ANOVA 
showed a marginally significant increase in V T 
(F = 4.363; P<0.05) after Zolpidem. There was some 
variation between individuals (see Figure 1), with 
Zolpidem causing a reduction in binding in one 
subject. 

[ 11 C]Flumazenil: For [ 11 C]flumazenil binding, a signi- 
ficant effect of Zolpidem was observed. Throughout 
the whole brain, V T was ~ 10% lower in the presence 
of Zolpidem (V T postplacebo: 4.10 ±0.56; postzolpi- 
dem: 3.70 ±0.37; F = 14.57; P< 0.0005; see Figure 2). 
The effect remained significant when all 83 regions 
were sampled (F = 123.7; P< 0.0001). No individual 
region survived correction for multiple comparisons. 
Zolpidem reduced [ 11 C]flumazenil V T in four of the 
six volunteers (see Figure 2B). 



Spectral Analysis: Assessment of Spectral Peaks 

The distribution of spectral peaks in different regions 
was identified using Clickfit. Each region yielded a 
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Figure 1 [ n C]Rol5-4513 scans. (A) Estimates of mean 
[ n C]Rol5-4513 Vj in different ROIs with and without Zolpidem. 
Error bars show s.e.m. of the five volunteers. (B) Intersubject 
variability, for the five volunteers, in [ n C]Rol5-4513 V T 
estimation in the whole brain. These data represent 
similar trends for each ROI within the subjects. ROI, region of 
interest. 



different spectrum of peaks, each peak position 
representing one of 100 values of fi, described in 
the 'Materials and methods' section. These results 
were pooled and a summed peak height for each 
peak was calculated for each region and each PET 
tracer in the presence or absence of Zolpidem. The 
average height for each /?-value was not representa- 
tive of the true peak distribution, because the peaks 
are discrete and the spectral analysis algorithms 
assign consecutive values for those peaks not 
matching a predefined 

There was a clear average peak height decrease 
with Zolpidem for both ligands. For ["C1R015-4513, 
the differences were seen between 0.0020 
and 0.0030 seconds 1 . With ["Qflumazenil, the 
decreases were seen at faster rates, between 0.0030 
and 0.0050 seconds 1 (data not shown). 

To further investigate the source of these differ- 
ences, a smaller number of brain regions with high 
and low al levels were selected for further analysis. 
Figures 3 and 4 illustrate data for the cerebellum and 
hippocampus after [ 1:1 C]Rol5-4513 and ["Cjflumaze- 
nil, respectively. As Zolpidem is a GABARal agonist, 
it was expected that regions with high al levels and 
low levels of other GABA receptor subtypes should 
show a greater change in peak profile than areas with 
lower al levels. In the cerebellum, a region with 
almost solely high al levels, the fast kinetic profile 
for [ 11 C]Rol5-4513 without Zolpidem represents al 
subtype binding (average height = 0.00306 seconds -1 ; 
range of peaks 0.0030 to 0.0040 seconds -1 , see Figure 
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Figure 2 [ n C]Flumazenil scans. (A) Estimates of mean fluma- 
zenil V T in different ROIs with and without Zolpidem. Error bars 
show s.e.m. of the six volunteers. In each region, there is an 
observable decrease in binding with Zolpidem administration, 
although significance did not survive post hoc correction. 
Significant trends were seen using a two-way ANOVA, with 
these combined regions, and all subregions in the whole brain. 
(B) Intersubject variability, for the six volunteers, in [ n C]fluma- 
zenil V T estimation in the whole brain. These data represent 
similar trends for every individual region, and more consistent 
decreases than for [ n C]Rol5-4513 lead to a more significant 
overall trend. ANOVA, analysis of variance; ROI, region of 
interest. 



3A). Zolpidem reduced average [ 11 C]Rol5-4513 peak 
height, lessening the contribution of those curves to 
the IRF and spreading their range (average 
height = 0.00185 seconds -1 ; range of peaks 0.0020 to 
0.0060 seconds -1 , Figure 3A). The effect of Zolpidem 
was estimated for this band of peaks in the 
cerebellum, to assess proposed al blockade. There 
was a significant decrease in Vs after Zolpidem 
(3.99±0.31 to 3.41 ±0.22, t=2.827; P<0.05). 

In contrast, in the hippocampus, a region with high 
oc5 levels but also containing the al subtype, there are 
two clear ranges of kinetic components of [ 11 C]Rol5- 
4513 (see Figure 3B). The fast component in the 
placebo scan (average height = 0.00122 seconds -1 ; 
range of peaks 0.0025 to 0.0045 seconds -1 ) is reduced 
and broadened in the presence of Zolpidem (average 
height = 0.00065 seconds -1 ; range of peaks 0.0020 to 
0.0060 seconds -1 ). This reduction in V s is significant 
(1.40 ±0.92 to 0.44 ±0.60, t= 2.369; P<0.05) and 
analogous to that found in the cerebellum, suggesting 
that this peak range reflected the al subtype. In the 
hippocampus, there was also a slower component, 
unaffected by Zolpidem (11.51 ± 3.44 to 10.00 ± 2.59, 
t= 0.8291; P> 0.45), which contributed to the majority 
of the convolution after Zolpidem administration 
(placebo: average height = 0.00179 seconds -1 ; range 



Journal of Cerebral Blood Flow & Metabolism (2012) 32, 731-744 



GABARzl contribution to [ 11 C]Rol5-4513 PET images 

JFM Myers era/ 



A 

0.006 
0.005 
0.004 
» 0.003 
0.002 
0.001 
0.000 



Cerebellum 



P r P r f 1 



Placebo 
Zolpidem 



737 



B 

0.006 
0.005 
0.004 
» 0.003 
0.002 
0.001 
0.000 



Hippocampus 



iiui 



-+JL 



Figure 3 Comparison of [ 11 C]Rol5-4513 spectral analysis results in (A) the cerebellum and (B) hippocampus, with data summed 
over five subjects. Only peaks 1 to 50 (frequency range 0.000631 to 0.00873 seconds -1 ) are shown. 
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Figure 4 Comparison of [ n C]flumazenil spectral analysis results in (A) the cerebellum and (B) hippocampus with data summed over 
five subjects. Only peaks 1 to 50 (frequency range 0.000631 to 0.00873 seconds -1 ) are shown. 



of peaks 0.0010 to 0.0014 seconds -1 ; Zolpidem: 
average height = 0.00177 seconds -1 ; range of peaks 
0.0010 to 0.0016 seconds -1 , Figure 3B). 

Comparable results for spectral analysis of [ 11 C]flu- 
mazenil in the cerebellum and hippocampus are 
shown in Figure 4. Although there are no clear 
distinctions between fast and slow kinetic compo- 
nents, there are differences in estimates of V T . For 



peaks between 10 and 50, there is a consistent and 
trend-level decrease in binding after Zolpidem 
administration in the cerebellum (6.11 ±1.18 to 
5.16 + 0.42, t= 3.343; P<0.06), and significant in 
the hippocampus (7.78 ±0.98 to 6.40 ±1.37, 
t= 2.923; P<0.05). The homogeneity of the spectra 
precluded any further separation of peaks based on 
regional receptor densities. 
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Occupancy Estimation 

The lack of effect on average V s by Zolpidem in the 
'slower' end of the spectrum, and the decrease in the 
'faster' range, led to broad estimates of i and j that 
would include all peaks in the two ranges. V s values 
were calculated for both slow range (i= 1; ;'= 25) and 
fast (i=25; 7 = 50), with and without Zolpidem. (V|) 
was plotted against (Vs— Vs) to produce occupancy 
plots, shown in Figure 5. For the slow range, little or 
no occupancy was shown, as expected because this 
was assumed to represent binding to the a5 subtype 
(Figure 5A). Two of the five regression plots were 
significantly different from zero, but with values of 
occupancy <20%, one positive and one negative. 
For the fast range, representing the al subtype, there 
are positive estimates of occupancy between 42.9 
and 84.8% (mean±s.d.: 64.3% ± 15.2%), with linear 
regression significantly different from 0 consistently 
with P< 0.0001 (Figure 5B). 

Occupancy was also estimated using plots for each 
of the six volunteers in the ["Clflumazenil cohort. 
The results are shown in Figure 6. V T was calculated 
by voxel-wise spectral analysis described earlier. 
Only 2 volunteers showed data fitted significantly by 
regression, giving Zolpidem occupancy estimates of 
30 and 33%. 



Discussion 

We have shown that Zolpidem significantly reduced 
["Clflumazenil V T in the whole brain by ~10%, 
whereas there is no reduction in [ 11 C]Rol5-4513 V T . 
These results suggest that the GABARal binding of 
[ 11 C]Rol5-4513 does not significantly affect overall 
V T , whether estimated by spectral analysis to pro- 
duce parametric maps of V T , or by compartmental 
modelling. However when frequency components of 
spectra derived from [ 11 C]Rol5-4513 images, pro- 
posed to reflect binding to the al subtype, are used to 
estimate V s , occupancy by Zolpidem is shown to 
range from 43% to 85%. The dose of Zolpidem was 
shown to be clinically relevant as it induced sedation 
and significantly impaired saccadic eye movement 
performance and word list recall. 

The difference in binding properties of nonselec- 
tive ligands in regions with differing populations 
of receptor subtypes ought to be demonstrable by 
applying various compartmental models. We hypo- 
thesised that more complex models would be 
required where there are multiple subtypes com- 
pared with brain regions with predominantly one 
subtype. One- and two-tissue compartmental models 
showed similar variability between individuals, 
but a different quality of fit to the data depending 
on region. In the cerebellum, an al-rich region, one- 
tissue compartment models provided the best fit for 
observed data, with higher and k 2 . Pharmacologi- 
cally, this corresponds to the rapid kinetic profile 
expected from the fast binding of a benzodiazepine 



ligand with relatively low affinity to the cerebellar 
subtypes. By contrast, in the hippocampus, the 
two-tissue compartment model was superior to the 
one-tissue compartment model. In addition, lower 
estimates of k 2 suggested a slower kinetic profile to 
that in the cerebellum. We suggest that this slower 
kinetic profile represents high-affinity binding to the 
GABARa5 subtype in a5-rich regions. These results 
are consistent with those of Maeda et al (2003) in 
nonhuman primates, who also showed that there are 
different kinetic components in regions containing 
different subtype levels and ratios. 

Both ["Clflumazenil and ["C]Rol5-4513 bind to 
more than one GABA a receptor subtype. Flumazenil 
binds with similar affinity (~lnmol/L) to al, a2, a3, 
and a5 GABA-benzodiazepine receptors and lower 
affinity at a4 and a6 subtypes (Sieghart, 1995). Rol5- 
4513 has 10- to 15-fold higher affinity for the <x5 
subtype (0.7nmol/L) compared with all remaining 
subtypes (Hadingham et al, 1993). Conventional 
compartmental modelling does not, in practical 
terms, accurately describe the complex pharmaco- 
dynamics when multiple receptor subtypes are 
present, especially when the different subtypes 
behave differently kinetically and are present in 
varying ratios depending on the brain region. In this 
study, compartmental modelling of ["C]Rol5-4513 
was not sensitive enough to detect changes in lower- 
affinity binding to al the subtype. Despite the 
usefulness of these models to show different regional 
kinetics, they are pharmacological simplifications 
and not suitable for accurately describing nonselec- 
tive ligands in quantitative analyses, although a 
two-tissue compartment model was proposed as 
effective in all regions by Asai et al (2009). Given 
that no one-compartmental model adequately repre- 
sented the kinetics of tracer binding in all tissues, 
and that there are multiple specific binding sites, we 
suggest this imposes limitations on using the 
reference region approach for a5-rich regions because 
it violates its underlying principle that the volume of 
distribution of nontarget bound ligand is the same in 
target and reference tissues (Lammertsma and Hume, 
1996). The pons is often used as a reference region, 
given the lack of suitable region devoid of benzodia- 
zepine receptors, but the low binding in the pons is 
not sufficient to fully meet the assumptions of the 
reference tissue model. Therefore, binding in regions 
of interest is underestimated and/or additional steps 
have to be taken to more accurately estimate such 
specific binding levels (Litton et al, 1994; Klumpers 
et al, 2008; Asai et al, 2009). 

In this study, spectral analysis consistently yielded 
higher estimates of V T relative to compartmental 
modelling. The whole-brain spectral profiles, using 
all 83 regions defined in the atlas (Appendix 1; 
Hammers et al, 2007), showed that with both 
["Clflumazenil and ["C]Rol5-4513, Zolpidem redu- 
ced the height of peaks in certain frequency ranges. 
Subtype differences could not be immediately 
inferred from these spectral frequency components, 
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Figure 5 Occupancy plots using the fractional volume of distribution in baseline and postdrug studies, according to equation (2), 
are shown for each subject. (A) Occupancy plots using l/ S i, defined by the 'slow' peak of the spectral analysis results (summing 
over 1 </c < 25). (B) Occupancy plots using l/ S 2, defined by the 'fast' peak (summing over 25</c<50). Zolpidem occupancy 
values obtained by linear regression, together with R 2 , are shown in each plot. Each data point represents the volume of 
distribution calculated for each ROI as listed in Appendix 1, excluding the ventricles, known to contain no GABA receptors. ROI, 
region of interest. 
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Figure 6 Occupancy plots for the [ 11 C]flumazenil cohort, showing V T calculated by voxel-wise spectral analysis in each region. 
Zolpidem occupancy values obtained by linear regression, together with R 2 are shown in each plot. Each data point represents the 
volume of distribution calculated for each ROI as listed in Appendix 1 , excluding the ventricles, known to contain no GABA receptors. 
ROI, region of interest. 



a problem compounded by complexity of the range 
of subtype ratios in the large number of brain 
regions included in the analysis. Unfortunately, the 
spectrum for [ 11 C]flumazenil proved too homo- 
geneous for separation of peaks based on regional 
receptor subtype density. However, for [ 11 C]Rol5- 
4513 images, in a5-rich regions such as the 
hippocampus and nucleus accumbens, we found a 
high-affinity binding site with a slower kinetic 
profile, consistent with the fitting of data from these 
regions with two-tissue compartment models. This 
slow component was not affected by Zolpidem 
administration, and was therefore attributed to 
high-affinity binding to the a5 subtype. In addition, 
in both a5-rich and a5-poor regions, there was a 
fast kinetic component present that was signifi- 
cantly reduced after Zolpidem administration, most 
likely representing specific binding with lower 
affinity to the al subtype. The intersubject variability 
of specific volume of distribution defined using 
these spectral bands was lower than for the total 
volume of distribution, V T . Therefore, spectral 



analysis was shown to be sensitive to detect the 
contribution of lower-affinity binding of [ 11 C]Rol5- 
4513 to GABARal in vivo. 

Confirmation of these subtype-specific frequency 
components in [ 1:1 C]Rol5-4513 images may be pro- 
vided with further study using a5-specific blockade. 
Until that is done, we cannot absolutely exclude that 
a2 and a3 subtypes may have contributed to the 
spectra we have taken to represent a5. In addition, 
there are no drugs selective for a2 and/or a3 
available. More broadly, these spectral band-pass 
methods may be useful to further sharpen analytical 
tools for in vivo PET assays, by reducing noise in the 
spectrum and by allowing possible subtype-specific 
probes into other nonselective ligands such as 
[ 11 C]diprenorphine and [ 11 C]PHNO. Similarly, sub- 
type-specific analyses of [ 11 C]Rol5-4513 PET in 
clinical and pharmacological applications may yield 
group effects hitherto overlooked because of the 
influence of other subtypes on the binding. 

Benzodiazepine agonists have previously been 
believed to exert their effects when occupying only 
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a relatively small number of receptors in humans. 
This is based on [ 11 C]flumazenil PET and [ 123 I]ioma- 
zenil SPET studies showing sedation or sleep is 
associated with clonazepam (0.03mg/kg) at 15% to 
23% occupancy, with diazepam (30 mg) at 24%, with 
alprazolam (2mg) at 16%, and with midazolam 
(50mcg/kg) at 17% to 35% or (6mg/h) at 20% to 
30% (Shinotoh et al, 1989; Pauli et al, 1991; 
Videbaek et al, 1993; Malizia et al, 1996; Fujita 
et al, 1999). Zolpidem is the only licensed benzo- 
diazepine subtype-selective drug with appropriate 
pharmacokinetics we were able to give to humans 
though, as it is an agonist, this limited the dose we 
were able to give safely. Abadie et al (1996) reported 
that 20 mg oral Zolpidem resulted in 20% to 29% 
occupancy of benzodiazepine receptors assessed 
with [ 11 C]flumazenil PET. We report that using 
[ 11 C]flumazenil, although occupancy cannot be 
directly estimated, the relatively small global 
decreases of ~10% in V T indicate this apparently 
low occupancy was made more indeterminate by the 
signal from receptor subtypes other than al. Lassen 
plots confirm that only two volunteers show signi- 
ficant regressions, estimating occupancy of ~30% 
with V NS >1.5. However, when using the higher 
frequency component of [ 11 C]Rol5-4513 to estimate 
V s at only the presumed al site, occupancy estimates 
ranged from 43% to 85%. This is accompanied by low 
V NS estimates, < 1 in most volunteers. What contri- 
butes to this range is unclear. We are unable 
to determine whether this was related to plasma 
levels because only samples from three indivi- 
duals were available and we only measured one time 
point and so we have limited pharmacokinetic 
information. Nevertheless, this higher occupancy is 
more consistent with that estimated by calculating 
that a 20 mg tablet given to a 100 kg individual will 
result in plasma levels of 650nmol/L. Assuming a 
volume of distribution in the brain of 0.5% and 90% 
protein binding, ~ 130 nmol/L will be available in the 
brain. As Zolpidem has an affinity of ~ 100 nmol/L at 
GABARal (Hanson et al, 2008; Hadingham et al, 
1993), the administered level of Zolpidem should 
result in high occupancy at this subtype. 

It was notable that for both PET ligands, there was 
high intersubject variability of the effect of Zolpidem 
on ligand binding. This is consistent with report by 
Abadie et al (1996) who also used 20 mg Zolpidem 
with ["Cjflumazenil PET and described 'spurious' 
intersubject variability including increases in 
[ 11 C]flumazenil binding in the presence of Zolpidem 
as we saw in some of our subjects. Therefore, it 
seems that subject heterogeneity is commonly seen 
when investigating the GABA system, of which one 
consequence has been some studies require large 
numbers of individuals. It is not clear what con- 
tributes to the variability of [ 11 C]flumazenil and 
["CjRolS^SlS binding but because both ligands 
bind to several different benzodiazepine subtypes, 
it is possible that relatively small variations in rece- 
ptor density ratios have a role. For instance, with 
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[ 11 C]Rol5-4513, the greater contribution to binding 
potential is proposed to come from the high-affinity 
GABARa5 component, which will be sensitive in 
ratio to the ubiquitous GABARal component. Simi- 
larly, it is possible that different subtypes are 
differentially affected by changes in endogenous 
GABA levels. It is for this reason that we have 
developed this technique with [ 11 C]Rol5-4513 to 
describe the contribution of benzodiazepine receptor 
subtypes to this variability. This method results in an 
importantly consistent spread of positive occupancy 
estimations across all individuals. 

We believe this development and application of 
spectral band-pass methods will help characterise 
benzodiazepine receptor subtypes as compounds 
with appropriate relative affinities become available. 
The use of nonselective ligands to show in vivo 
occupancy of receptor-selective compounds will be 
more accurate as specific binding of the ligand is 
estimated. At present, [ 11 C]flumazenil PET has been 
used to determine occupancy by two such com- 
pounds, TPA023 at a2 and a3 and a5lA at a5 (Atack 
et al, 2011). As [ 11 C]flumazenil is nonspecific and 
a5lA and TPA023 each have equivalent affinity for 
al, a2, a3, and a5 subtypes, the reduction in 
[ 11 C]flumazenil was presumed to reflect occupancy 
equally at all subtypes to infer occupancy at a2/a3 for 
TPA023 and a5 for a5lA. Eng et al (2010) emphasised 
the importance of establishing the plasma-occupancy 
relationship in humans rather than merely relying on 
preclinical data because in both studies, the EC 50 in 
humans was lower compared with preclinical data. 
Therefore, the approach we have adopted here will 
greatly enhance our ability to understand the role of 
different benzodiazepine subtypes and its potential 
to determine contribution of receptor subtypes in 
other nonspecific ligands. 
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Appendix 1 
Region List 



(R = right, L = left) 
R hippocampus 
L hippocampus 
R amygdala 
L amygdala 

R anterior temporal lohe — medial 

L anterior temporal lobe — medial 

R anterior temporal lobe — lateral 

L anterior temporal lobe — lateral 

R parahippocampal and ambient gyri 

L parahippocampal and ambient gyri 

R superior temporal gyrus 

L superior temporal gyrus 

R middle and inferior temporal gyri 

L middle and inferior temporal gyri 

R lateral occipitotemporal gyrus (fusiform gyrus) 

L lateral occipitotemporal gyrus (fusiform gyrus) 

R cerebellum 

L cerebellum 

Brainstem 

L insula 

R insula 

L occipital lobe — lateral 

R occipital lobe — lateral 

L gyrus cinguli — anterior part 

R gyrus cinguli — anterior part 

L gyrus cinguli — posterior part 

R gyrus cinguli — posterior part 

L frontal lobe 

R frontal lobe 

L posterior temporal lobe 

R posterior temporal lobe 

L parietal lobe 



R parietal lobe 
L caudate nucleus 
R caudate nucleus 
L nucleus accumbens 
R nucleus accumbens 
L putamen 
R putamen 
L thalamus 
R thalamus 
L pallidum 
R pallidum 
Corpus callosum 

R lateral ventricle — frontal horn, central part and occipital 
horn 

L lateral ventricle — frontal horn, central part and occipital horn 
R lateral ventricle — temporal horn 
L lateral ventricle — temporal horn 
Third ventricle 

L precenntral gyrus — frontal lobe 
R precenntral gyrus — frontal lobe 
L straight gyrus — frontal lobe 
R straight gyrus — frontal lobe 
L anterior orbital gyrus — frontal lobe 
R anterior orbital gyrus — frontal lobe 
L inferior frontal gyrus — frontal lobe 
R inferior frontal gyrus — frontal lobe 
L superior frontal gyrus — frontal lobe 
R superior frontal gyrus — frontal lobe 
L postcentral gyrus — parietal lobe 
R postcentral gyrus — parietal lobe 
L superior parietal gyrus — parietal lobe 
R superior parietal gyrus — parietal lobe 
L lingual gyrus — occipital lobe 
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R lingual gyrus — occipital lobe 
L cuneus — occipital lobe 
R cuneus — occipital lobe 
L medial orbital gyrus — frontal lobe 
R medial orbital gyms — frontal lobe 
L lateral orbital gyms — frontal lobe 
R lateral orbital gyrus — frontal lobe 
L posterior orbital gyrus — frontal lobe 
R posterior orbital gyrus — frontal lobe 
L substantia nigra 



R substantia nigra 

L subgenual frontal cortex 

R subgenual frontal cortex 

L subcallosal area 

R subcallosal area 

L presubgenual frontal cortex 

R presubgenual frontal cortex 

L superior temporal gyrus, anterior part 

R superior temporal gyrus, anterior part 
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